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The past decade has witnessed a great deal of work on the theoretical 
modeling of spin-glasses. The efforts have been mainly concentrated on the 
study of the equilibrium properties of mean-field and low- dimensional models. 

Thanks to the the replica method there is at present a good understanding 
of the Gibbs-Boltzmann measure in mean- field models JTJ. The phase-space 
turns out to have a very complicated landscape with many hierarchically 
organised states. These states are separated by barriers diverging with the 
size of the system and hence implying ergodicity breaking. 

The implementation of the replica approach for low- dimensional models 
has proven to be a hard problem and few analytical results are at present 
available for them. Whether the picture found from the replica analysis of 
mean-field models can be applied to realistic systems is still an open question; 
the nature of their ordered state and phase transition remain controversial 0, 
|3|] . Indeed, it is sometimes argued that mean-field models are too unrealistic 
and cannot be accepted as candidates to describe real spin glasses. 

However, the experimental evidence suggests that spin-glass physics in 
the low temperature phase is essentially dynamical: they continue to evolve 
long after thermalization at a subcritical temperature and their dynamical 
properties depend on their history throughout the experimental time. These 
effects are known as aging phenomena and are widely observed in many other 
disordered systems. Evidently, the relevance of a model consists in its ability 
to capture these phenomena. 

Besides, the experiments suggest that the dynamics of spin-glasses has 
similar features to those found in usual glasses ||. One may then hope to 
apply the developments in the (perhaps simpler) area of spin-glasses to other 
glassy systems. 

Though several phenomenological approaches to aging phenomena in spin- 
glasses have been made 0, a fully microscopic description is still lacking. It 
was recently noticed that mean-field models, though schematic, do possess 
a rich dynamical phenomenology resembling, at least qualitatively, the ex- 
perimental observations 0. Furthermore, it turns out that one can obtain 
analytical results for their long-time dynamics. 

In what follows we shall describe the long-time dynamical behaviour of 
two representative examples: the p-spin spherical model defined by 





(1) 
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and the Sherrington-Kirkpatrick (soft-spin) model 

N / N \ 

Hj[s] = -T, J^ + u £(^-l) 2 • (2) 

The last terms enforce the spin measure (spherically constrained and soft 
spin, respectively) , and the quenched disorder J has in both cases zero mean 
and variance J =p\/(AN^) and J = 1/(2JV), respectively 

The analytical study of the long-time relaxation of these models reveals 
that both models have aging phenomena, but with rather different char- 
acteristics 0, Their behaviour can be interpreted roughly on the lines 
introduced by Bouchaud ||. 



1 Geometry of phase-space 

As a first approach to the dynamical problem it is convenient to have an 
understanding of the low-temperature structure of phase-space. We have 
at hand two tools: the replica solution and the Thouless-Anderson-Palmer 
(TAP) approach. The former gives us very detailed information of the lowest- 
lying ('equilibrium') states which contribute to the Gibbs-Boltzmann mea- 
sure, while the latter gives us an idea of what the landscape looks like for 
higher free-energies. Both approaches reveal that the two models we are con- 
sidering have some important qualitative differences, as the dynamics also 
shows. 

In the case of the p-spin spherical model, one level of replica symmetry 
breaking is exact. The equilibrium states are separated by O(N) barriers 
(the difference in energy and free-energy between them being 0(1)) []TT|] . 
The 'size' of the equilibrium states is q^A while the overlap between different 
states is q = |12| . 

The TAP equations for this model have solutions for a free-energy range 
(F lrsb , F thres ), corresponding to an energy range (E lrsb , E thres ). The thresh- 
old energy E thres (free-energy F thres ) is greater than the equilibrium energy 
E lrsb (free-energy F lrsb ) by a (temperature and p dependent) 0(N) value. 
Under the threshold the local minima of the TAP free-energy are separated 
by 0(N) barriers while above the threshold there are no minima at all [|ll| 



The free-energy Hessian of a solution has eigenvalues larger than a certain 
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Xmin that depends on the energy of the solution. For sub-threshold solutions 
A m in is positive-definite, its value decreases with increasing free-energy until 
it vanishes at the threshold. As the temperature is changed the solutions 
neither merge nor bifurcate, and their free energy changes smoothly. 

In the case of the SK model, the equilibrium states are also organised in 
an ultrametric fashion, but with all possible distances between them allowed. 
The barriers between equilibrium states are believed to be of only 0(N a ) with 
a ~ 1/3. The solutions of the TAP equations tend to split as the temperature 
is lowered in a second-order fashion. The Hessian of the solutions has a 



spectrum going down to zero |13[], so that there is the possibility that the 
barriers between neighbouring high free-energy TAP solutions be finite. For 
this model there seems to be no high free-energy threshold below which all 
TAP solutions are separated by infinite barriers and above which there are 
no solutions, as happens in the p-spin spherical model. 

Fig. 1 shows schematically this structure of states and barriers. 



2 Asymptotic out of equilibrium regime 

Let us now turn to the dynamics of an infinite system. We take the ther- 
modynamic limit N — ► oo from the outset [], large times henceforth mean 
t —>■ oo after N — > oo. Under this hypothesis there may be unsurmountable 
barriers. We consider a process in which the system evolves from an initial 
configuration chosen at random. 

Taking into account the structure of the phase-space, three possibilities 
for the long-time dynamics can be expected: 

A. The system reaches the Gibbs-Boltzmann distribution in a character- 
istic time t eq . 

B. The system equilibrates within a separate ergodic component. The 
distribution is a Gibbs-Boltzmann distribution restricted to that sector of 
phase-space, and it does not evolve after a characteristic time t eg . 

C. The system does not reach a time-independent Gibbs-Boltzmann dis- 
tribution within any fixed sector of phase-space at any time. There is no t eq 

1 This is the main difference between this dynamicai approach and the one proposed 
by Sompohnsky Q. 
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such that for t > t eq the distribution has stabilised. 

Case A corresponds to an ordinary equilibration process, thermodynami- 
cal and long-time dynamical calculations coincide. The equilibrium theorems 
hold for times larger than t eq . 

Case B corresponds to falling in a high-lying stable state and a purely 
statical calculation does not necessarily yield the right values for order pa- 
rameters and free-energy. However, from the dynamical point of view this 
situation is not very different from the situation A: after the time t eq the 
system is for all practical (dynamical) purposes as if in equilibrium. Time- 
homogeneity and the fluctuation-dissipation theorems (FDT) hold for any 
two times t\ > t 2 > t eq . 

Case C is what we find in mean-field spin-glass models. As the system 
evolves, the dynamical free-energy density decreases, and the system finds 
an ever decreasing portion of the landscape where it can move. While at 
high temperatures this process ultimately leads the system to the equilibrium 
state, at low temperatures the geometry of the region of phase-space available 
at a certain time becomes more and more complicated, and the system is 
slowed down and does not reach an equilibrium state in finite times. However, 
it does not remain permanently confined in any finite region within which it 
is in local equilibrium. 

For instance, in the scenario of Bouchaud ||, one considers a phase- 
space with 'traps' separated by finite barriers with a wide distribution of 
lifetimes. Traps visited at long times tend to have long lifetimes, for purely 
probabilistic reasons. The system does not reach the true states (those with 
infinite lifetime) in finite times, hence the name ir weak ergodicity breaking'. 

Interestingly enough, once the thermodynamic limit has been taken there 
are no external parameters in the dynamics of the system that go to infinity. 
Instead it is the age t w (the time elapsed since the quench from above the 
transition temperature) that regulates the slowing down of the system. For 
longer waiting times, the system has the opportunity of finding deeper 'traps' 
and narrower channels, it becomes less susceptible to changes and thus it 
ages. 
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3 Assumptions on the Correlation and Re- 
sponse Functions 

The solution of the models is based on some assumptions that take into 
account the preceding scenario. We now describe these assumptions. 

i. 'Weak-ergodicity breaking' 

After any time t w (that may be interpreted as a waiting time) the system 
continues to drift away and it reaches, asymptotically, the maximum distance 
compatible with the remanent magnetization. Thus, the correlation function 
C(t + t w: t w ) = (1/N) E?=i (si(t + t w )si(t w )) satisfies 

dC(t + t w ,t w ) dC(t,t') 

dt - u - u {i) 

t > t', and in the absence of a magnetic field 

lim C(t + t w , = V fixed t w (4) 

t—>oo 

This hypothesis is supported by the numerical simulations of various 
mean-field || [T7|, [18| and realistic |19| models. The numerical solution of 
the mean-field dynamical equations for the p-spin spherical model (see Ref. 
[0]) also support this assumption: in Fig. 2. a. we plot the decay of the 
correlation function C(t + t w ,t w ) vs. t + t w . It is clear that the sign of the 
derivative is negative for the whole time-interval plotted. Moreover, Fig. 2.b. 
shows the decay of the correlation C(t, £') as a funtion of t' and the derivative 
is positive in this case. (T = .3, p = 3 in both cases.) 

ii. 'Weak long-term memory' 

The response to a constant (small) magnetic field applied from t' = up 
to if = t w , h tw (t') = h9(t w — t')9(t'), i.e. the thermoremanent magnetization, 
decays to zero for a long-enough time t 

rt w 

lim m tw (t) = lim / dt'G(t,t') = V fixed t w (5) 

t— >00 ™ t—rOO JQ 



2 We present these figures to show the qualitative tendency of the correlation and re- 
sponse functions, a more precise numerical analysis is in order to obtain quantitative 
results. 
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The associated experimental curves [^0], and the simulations of the D- 
dimensional spin-glass cell are compatible with this assumption. 

It is this property that makes solutions to the long-time dynamics possi- 
ble: Because the memory of the system of its history is strong when integrated 
over long times but weak at each given time, one does not need to solve all 
the short time (far off-equilibrium) details. 

in. Fast and slow relaxations. 

After a (long) time t' there is a quick relaxation in a further time r = t — t' 
(r short compared to t') to some value q (the dynamical Edwards- Anderson 
parameter), followed by a slower drift away. Fast and slow relaxations cor- 
respond to relaxations within and away from a trap, respectively. Within 
a trap the system behaves as if it were in a local equilibrium, the equilib- 
rium properties (Fluctuation-Dissipation Theorem (FDT) and homogeneity 
in time) are assumed to hold. The value q is the size of the traps or the 
width of the channels through which the system evolves. 



C(t,t') = C FDT (t-t')+C(t,t>) 

G(t,t') = G FDT (t-t') + g(t,t') 



GFDT\t — t') 



dC FDT (t-t') 
dt' 



with 



CW(0) = l-q lim t _^ 0O CW(t-t') = 
C{t,t) = q lmw^M') = 0. 

The relaxation away from a trap is slow: 

dC(t,t') dC(t,t') 

— -A-^~0 for large t,t' (6) 

In the simulations of the D-dimensional hypercubic cell of Ref. ||, the 
log(C(t + t w , t w )) vs. log(t) curves show that for short times t = (t + t w ) —t w 
compared to t + t w all the curves for different waiting times t w lie on each 
other. Hence, the correlation function is homogeneous in time for that range. 
The fast relaxation to a value q can be checked in the curves presented in 
Ref. 0. 

In the numerical solution of the mean-field dynamical equations of the 
p-spin spherical model we also find support for this hypothesis. The Figs. 
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2. a and 2.b show a 'fast' decay from the value 1 at equal times to a certain 
value. This is more clear in Fig. 2.b. where C(t,t') decays fastly from 1 to 
a value q ~ .9 This is justified since the total time involved is larger and the 
interval for the FDT decay must be short compared to the total time but it 
itself must be long. 

Two further properties that are in part consequences of the assumptions 
already made are: 

iv. For long times t,t' C and Q are related by 



where X depends on the times only through C. 

The relation (|7|) with X[t, t'\ = X[t,C(t,t')] implies no assumption. One 
can show that in the large times limit the explicit dependence on t can be 
neglected if the energy, the susceptibility, and all the 'generalised suscepti- 
bilities' || 0] have a limit. 

If we supplement the definition of X[z] with X[z\ = 1 for q < z < 1 then 
the relation (|7|) holds for all C(t,t'), t' large. X[z] may be discontinuous in 
z = q, where it jumps from X[q] to 1 (X = 1 ^ FDT). 

An assumption concerning the function X(C), suggested by the simula- 
tions is: 

v. X[C] is an increasing function of C. 

vi. 'Triangle relations' between the correlations at three large times i TOl - n < 

tint ^ t max 

The monotonicity of the correlation function with respect to both times 
described above, allows to write 



We can consider the limit for large t min , and for fixed C{t maxi t min ). The 
fact that such a limit exists implies that we can write, for large enough t m i n : 



g(t,t') = x[c(t,t')} 



dC(t,t') 
dt' 



9(t - 1') 



(7) 




(8) 




(9) 
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where the explicit dependence on times have dissapeared, and the three cor- 
relations are related by 'triangle' relations. It is convenient to define the 
inverse relation: 



C '(t max ; tint) f \C(t max j^min)) C(^min) ^ini)] (10) 

A simple computation involving four times shows that the function / is 
associative [0] 

/(/(gi,? 2 ),?3) = /(?i,/(?2,g3)) (ii) 

The property of weak ergodicity breaking implies that 

f(Qi,Q2) < mm(q 1 ,q 2 ) (12) 

One can now classify all the possible associative / that satisfy ([12]) as follows: 
consider the values q* for which 

MA*)=<fi- (13) 

We call them fixed points. Fixed points and intervals between consecutive 
fixed points are 'correlation scales'. The latter we call 'discrete scales' and, 
in particular, the interval (q, 1) is one of them. A dense set of q* corresponds 
to a dense set of scales. 

It is shown in Ref. [0] that if q\ and q 2 belong to different scales then 

fill, <h) = min(<7i, q 2 ) (14) 
If, instead, they belong to the same discrete scale: 

f(Qi, 92) = 3 k 1 bk(qi)jk(q2)] , (15) 

where jk are functions that may be different for each discrete scale. 

In Fig. 3 we sketch the function f(a, a) for a ferromagnet, a paramagnet, 
the p-spin spherical spin-glass and the SK spin-glass. In the next sections 
we present the analysis of the p-spin spherical model and we discuss the 
corresponding results for ths SK model 0. 
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4 A Simple Example 



The analysis and assumptions in the preceding section are quite general and 
model-independent. The precise form of X(C), the limits of the correlation 
scales, and the functions jk associated with the discrete scales (if any) have 
to be extracted from the dynamical equations of each model. 

The analysis of the SK model along these lines was done in Ref. [0], 
here we work out the dynamics of the p-spin spherical model within this 
formalism. We shall afterwards discuss the differences between these two 
representative models. 

The mean-field equations of motion for the p-spin spherical model read: 

dC(t,t') 



Of 



dG{t, t') 



- (1 - p/3S(t)) C(t, if) + 2 G(t', t) 

+ fi f dt"C p - 1 (t,t")G(t',t") 
Jo 

+ fi(p-l) t dt" G(t, t") C p ~ 2 (t, t") C(t", t') 
Jo 

-{l-p(3E{t))G{t,if) + 5{t-if) 

+ (i(p-l) f dt" G(t, t") C p ~ 2 (t, t") G{t", if) 
Jf 



where /i = p{3 2 /2 and the energy E: 



E(t) 



dt"C p -\t,t")G(t,t") 



(16) 



(17) 



(18) 



Considering assumption Hi we can separate the small time difference 
(FDT) regime from the regime of widely separated times. We hence obtain: 



+ 1)C fdt (t) + (v+pP£oc) (l-C FDT (r)) =fi dr" C^r-r" 



dC- 



dr" 



fdt V m 



with the asymptotic energy Eoo given by 
_P r 
2 



i- q p )+p f t dt"g{t,t")c p ~ 1 (t,f) 

Jo 



The correlation decays to a value q determined by 
-(l-p(3E OQ )+n(l-q p - 1 ) = - 



(19) 
(20) 

(21) 
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Equation ([19] 

the relaxation 'within a trap 



corresponds to the dynamics a la Sompolinsky-Zippelius of 



For the long time differences we obtain: 







g(t, t>) [-(1 - q)- 1 + fi (1 - q)(p - 1) C?- 2 (t, t') 

+ti ( P - 1) f dt" g(t, t") c p - 2 (t, t") g(t", t') , 

Jo 

C(t,t') [_(l-g)-l +,x(l- g ) C^^t') 

+fi f dt"c p -\t,t")g(t',t") 

Jo 

+0 (p - 1) /" dt" g(t, t") c p - 2 (t, t") c(t", f) . 

Jo 



(22) 



(23) 



We have now two (coupled) sets of equations, corresponding to two time 
regimes. It is important to note that in order to claim that eqs. (|2"2"D and 



( |23|) are asymptotically valid we need the property of weakness of the long- 
term memory, which allows us to disregard in the integrals any finite-time 
interval. We can then neglect the fact that for the initial times the asymptotic 
equations do not hold. 

Above the critical temperature g ~ 0, the long-term memory is absent, 
and we are left with only equation ([19]) [15] . This is the same equation that 



arises in the mode-coupling theory of structural glasses |Id[ . 

In the glass phase, in order to determine any quantity, we have to solve 
the whole set of eqs. (JT|), (|22|) and (g|). 

We concentrate here on the solution of eqs. 



(H) and (H). We have 
neglected in them the time derivatives following in. This brings as a conse- 
quence the fact that these equations (but not eq. (|19D) have a continuous 
set of reparametrization invariances. Indeed, from a solution C, g we obtain 
infinitely many others: 



c(M') = c(h(t),h(t')) , g(i,i r ) = h'{t')g{h{t),h{t')) 



(24) 



with h any increasing function. This invariance is a consequence of hav- 
ing neglected the time derivatives in making the asymptotic limit. The full 
dynamical equations have no such invariances; because of causality their solu- 
tion is unique. We face the selection problem, a rather common phenomenon 
in the asymptotic limit of solutions of differential equations. In this work we 



10 



discuss solutions only modulo reparametrizations, the selection problem for 
this kind of dynamics remains to be solved. 

Let us now use the assumptions iv and vi to write ( |2"2"D and fl2"3| ) in terms 
of X and /. Defining 

F[C] = - £dC'X[C'} H[C} = - jT dC'C' p ' 2 X[C] (25) 

the dynamical equations ( P2"| ) and ( |21f ) become 
= -{l-q)^F[C]+pL{p-l){l-q)H[C\ 

+/*(p- 1) /Vx[C']C'^ 2 F[7(C',C)] , (26) 



.(1 _ g )-i+ A* (l-?)d>- a 



C-/i(p-l) / dC'C*- 2 F|7(C,C')] 



£(f) 



+//(p-i) / rfc'x[C]c^V(c,c') 

JO 

+/i(p-i) ^dCx[C]c*- 2 7{C,c), 

'l-q p )+p I" dC XIC'IC'?- 1 
Jo 



2 



(27) 
(28) 



In using the relation (0) and (§) we have eliminated the times. Hence 
we have divided by the reparametrization group. Indeed, (|7|) and (|9|) are 
reparametrization-invariant relations themselves. 

The next step is to investigate the possibility of having 'discrete' correla- 
tion scales. The method is to propose that there is a scale, and then check 
which limits the equations allow the scale to have. Using the properties of 
the function / described in the previous section 0, these equations can be 
rewritten within a discrete scale C G (a^, a*), < a* < q, 







1 _ F[C ] + p ( p _ 1) (1 _ q ) H [C] - /x (p - 1) F\C] H[a\] 



C 



dC'X[C'}C' p - 2 F[f{C',C)} , 
l- q )-i +fI (l- q )C p ~ 2 -v(p-l)H[al} 

fi(p-l) r dCC >p - 2 (C'X[C] - F[C'\) 
Jo 

[dC'C'P- 2 (X[C')7(C,C')-F[J(C,C')]) 



(29) 



+ 

+fl (p - 1 
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+A*(p-1) f al dC'X[C')C'P- 2 f{C',C). (30) 

Evaluating eq.(|2D) in C = a\ we have 

- (1 - q)- l F[a\} + //(p - 1) ((1 - g) - F[a*]) F[a*] = , (31) 
and differentiating eq.(^) w.r.t. C and evaluating in C = a\ we get 

(1 - q) (<- 2 - q?- 2 ) - a^ 2 F[al] - H[a{] = , (32) 
provided X[a*] ^ 0. From eqs. (|3l|) and (|3^) F[a*] is given by 



J? M = ( 1 -9)( 1 -y^l) • ( 33 ) 

Differentiating again eq. (|29|) w.r.t. C, evaluating the result in a* and using 
the previous eqs. (j32|) and ([33|) we finally get 



Since we assume that X[C] is an increasing function of C (assumption v) then 

a* = g (35) 

and 

XM=X[ g ]= (p " 2)(1 " g) . (36) 

Once we have shown that the upper limit of the discrete scale must be 
q we need to obtain its lower limit. First, we shall show that X[C] must be 
constant inside the scale, i.e. X[C] = X = (p — 2)(1 — q)/q, VC G (a^q). 
Computing the derivatives of eqs. ( p9|) and ( |30"D w.r.t. C, multiplying the 
former by X[C] and subtracting we have 



= f° dC C' p - 2 X[C] fiC 'P (x[f(C,C)] -X[C'} 

+ £ dC' C /p - 2 X[C] 2^ (X[7(C, C)\ - X[C}) . (37) 
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This equation has X[C] = X as a solution. To show that it is the only 
admisible solution we derivate again w.r.t. C and evaluate in C = a* = q to 
get 

^1 = =► X[C]=X, (38) 

VCG (a* 2 ,q). 

Finally, using the constancy of X we calculate F[C] and H[C] 

F[C] = -X(q - C) H[C] = ~^-{q p - 1 - , (39) 

p — 1 



and from eq. fl2"5|) we get 

a ; = or a; = g . (40) 

Hence the discrete scale is 'empty', there are no discrete scales at all but 
a continuos set of fixed points from to q, or there only one discrete scale 
that spans the whole interval (0, g) (apart from the FDT sector). We have 
checked that the former possibility is excluded by the dynamical equations, 
namely (f27|). 

We can now turn to the explicit solution within the discrete scale spanning 
the interval (a^ = 0, a* = q). Eqs. (p9[) and ( p0|) can now be written as 

= g - C - p (p - 1) (1 - g) 2 fa*" 1 - ff^ 1 ) 

+^(p-l)(l-g)X l c 9 dC'C' p ~ 2 J(C',C), (41) 

= C[-l+/i(p-l)(l-g) 2 C p - 2 " 

+/x (p - 1) (1 - g) X f c dC C' p - 2 J(C, C) . (42) 

The value g can be immediately determined from these equations; they imply 

1 



p — 1 



W p ~ 2 (l ~ qf ■ (43) 



From the analysis presented in Ref. 0] we know that the solution within a 
discrete scale has the form 



(44) 
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C < C . We should now use the equations ([41]) and ( f42|) to determine the 
function J(C',C). 

Making the change of variables: 

P = ln (fj ! Kp)=3-\qe p ) (45) 
eqs. (|41|) and (|2[) both become [] 

[1 - F~ 2 (p)] + ^ T dp' 9 Jf^l k{p - p ')=0. (46) 
p — 1 jo op 

The smooth solutions to this equation are: 

k(p)=e^ =► 7( C ',C) = g- 1 ^ (47) 
Which in turn implies that C is of the form: 

c<«,0=,$J («) 

for any increasing /i. 27ws is as /ar as w;e can go analytically without solving 
the selection problem. The numerical solutions suggest || that h(t) = t 1 . 
For large t w one has, as t — t w — > oo 

C(t + ^,tJ^ (M^P) (49) 



which explicitly shows aging. 



5 Discussion 

We have shown how to obtain some analytical results for spin-glass mean-field 
dynamics, on the basis of the assumptions made. In this solution the aging 
phenomena are explicit. The solution is not complete, some quantities we 
have only obtained modulo time-reparametrizations. However, the long-time 

3 Remarkably, this is the same equation that one obtains for long times in the Som- 
polinsky (time-homogeneous) dynamics Jig] with p playing the role of time-differences. 
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limits of all quantities that depend on a single time (energy, magnetization, 
q, etc) are unaffected by reparametrizations, and hence fully determined. 

We have found a solution for the p-spin spherical model with one discrete 
scale apart from the FDT scale. Let us now briefly describe what happens 
when one applies exactly the same procedure to the SK model ]7|. It is 
found there that the dynamical equations do not admit any discrete scales 
(apart from the FDT (q, 1) scale), but only a dense set of scales in the 
correlation interval (0,g). This is very much like in the static treatment, 
discrete scales playing the same role as levels of replica symmetry breaking. 
There is however a difference: while for the SK model one finds that the 
values of the asymptotic energy, q, and transition temperature coincide (to 
O(N)) with those obtained in the static treatment, for the p-spin model 
they do not. Indeed, for this model one finds that the values for q and E 
correspond to those of the threshold level. 

Another difference between the two models is in the behaviour of the 
thermoremanent magnetization. In the p-spin spherical model it decays with 
a rate that is inversely proportional to the waiting time the field has been on. 
In the SK model the decay of the magnetization after a long waiting time 
can be seen as taking place in steps, each of which is takes much longer than 
the preceding one. Hence, the SK model has a larger degree of 'freezing' of 
the magnetization. 

To conclude, let us remark that the dynamical spin-glass phase can be 
viewed as a phase where a symmetry is spontaneously broken. The dynamical 
process has a supersymmetry (SUSY) group of invariances [22| with three 
generators 

D' ; IF ; P? ,D 1+ = J- + J- (50) 

associated with probability conservation, FDT and time-homogeneity. The 
two-time (t, t') functions are defined in the time region whose two boundaries 
are the lines t = and t' = 0, respectively. 

The full SUSY group is broken down in certain time-regions down to the 
subgroup generated by only D'. In the high temperature phase the effect of 
the initial (boundary) condition is to break the SUSY in a region of width 
~ t eq near the boundaries. In the low temperature phase the SUSY-breaking 
persists in an infinite region away from the boundaries. This is much like the 
effect of symmetry-breaking boundary conditions in ordinary sytems. 
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